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Linear perturbations of nomothetic self-similar stiff fluid solutions, S[n], with circular symmetry 
in 2 + f gravity are studied. It is found that, except for those with n = 1 and n = 3, none of them 
is stable and all have more than one unstable mode. Hence, none of these solutions can be critical. 
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O . I. INTRODUCTION 
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\Q • Studies of the non-linearity of the Einstein field equations near the threshold of black hole formation reveal very 
rich phenomena [1], which are quite similar to critical phenomena in statistical mechanics and quantum field theory 
[2]. These phenomena were first discovered by Choptuik in a study of a collapsing massless scalar field, and since that 
■ time they have been found to exist in various other matter fields [3]. 

In this paper, we study the gravitational collapse of a circularly symmetric stiff fluid with homothetic self-similarity 
(HSS) in 2 + 1 gravity. It is well known that a time-like massless scalar field is energetically equivalent to a stiff fluid 
| [6]. In fact, introducing the quantities 
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where <f> denotes the scalar field and ( ) tfl = d( )/dx fi . It should be emphasized that the above equivalence is valid 
only when the massless scalar field is time-like. 

The main motivation for the present study comes from a recent investigation of the critical collapse of a stiff fluid in 
3+1 gravity [7] , in which it was shown numerically that a massless scalar field and a stiff fluid exhibit quite different 
critical phenomena near the threshold of black hole formation, given the equivalence of Eq. (1.1). In particular, the 
critical solution of the scalar field has discrete self-similarity (DSS), and the corresponding exponent 7 is 7 ~ 0.374 
[1], while for a stiff fluid, the critical solution has HSS, and the corresponding exponent is 7 ~ 0.94 [8]. It should be 
noted that in the latter case, the HSS (critical) solution is not stable with respect to kink perturbations [9]. 

In this paper, we study the same problem, but in 2 + 1 gravity. In particular, we first review some properties of 
the self-similar solutions in terms of a stiff fluid particularly focusing on their extension across the self-similar horizon 
(SSH). This is crucial when we consider the boundary conditions for their linear perturbations. As we show, in some 
cases the extension for the scalar field is analytical, and the only change is that the scalar field becomes space-like 
in the extended region, while for a stiff fluid such an extension leads to a negative energy density and space-like 
velocity. Thus, for the latter, we have to drop the requirement of analytical extension. Because the background is not 
analytical, there is no reason to require the perturbations to be analytical across the SSH. As we show, removing the 
requirement of analyticity results in a different spectrum for the perturbations. In particular, we find that, except 
for the solutions with n = 1 and n = 3, all the self-similar solutions for a stiff fluid are unstable and have more than 
one unstable mode. Hence, none of the self-similar solutions is critical for the gravitational collapse of a stiff fluid, 
because, by definition, a critical solution has one and only one unstable mode. 

It should be noted that, because the problem has been studied extensively by several authors, it is difficult to aviod 
repeating material that has previously been presented, but we try to keep such repetition to a minimum. In Ref. 
4) and 5), only linear perturbations of the self-similar solutions that have analytical extensions across the SSH were 
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studied, although some of those results are valid even for solutions that do not have analytical extension. In this 
paper, we study perturbations of all self-similar solutions, both those with and those without analytical extension, in 
terms of a stiff fluid. 



II. SPACETIMES WITH HOMOTHETIC SELF-SIMILARITY 



The self-similar solutions of a stiff fluid are given by [5,10] 
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where po = 2 4c c 2 and n = 1/[2(1 — c 2 )], with n (or c) being a free parameter. The above solutions in general are 
valid only in the region satisfying u < 0, v < and v > u in the (u, w)-plane, which is referred to as Region II and 
is shown in Fig. 1. Region /, in which u < 0, v > 0, \u\ > v, is considered an extended region. The extension is 
analytical only when n is an integer. However, when n is an even number, Eq. (2.2) shows that the energy density p 
becomes negative in Region I. This is not physically feasible. Given this situation, one may be tempted to consider 
the extension obtained by replacing — v with \v\. However, this extension is not analytical. In fact, no analytical 
and physically reasonable extension exists for the stiff fluid solutions of Eq. (2.2), except in the case that n is an odd 
integer. This fact is crucial when we consider the boundary conditions at v — of linear perturbations. 

It should be noted that the surface v = is a sonic line of the stiff fluid. This surface is also called a self-similar 
horizon, that is, a null surface of constant z, defined in Eq. (3.2) [7]. In addition, this SSH is degenerate, as shown 
explicitly in Refs. 5) and 10). This differs from the corresponding (3+ l)-dimensional case [7], in which the apparent 
horizon forms after the SSH. For this reason, the boundary conditions for the linear perturbations in these two cases 
are different. This, in turn, affects the spectra of the perturbations. 



FIG. 1. The Penrose diagram for the solutions given by Eq. (2.2). Region I is the extended region. When n = 2rn + 1, 
where m is an integer, the dashed line v = (or v = 0) represents an apparent horizon, and the spacetime is singular on the 
horizontal double line r = (— u) 2n — (— v) 2n — 0. For n / 2m + 1, the extension is not analytical. 



Let us first note that, when we consider perturbations of the self-similar solutions of Eq. (2.2), we consider only 
Region i7, where the equivalence (1.1) between the stiff fluid and the massless scalar field is always valid, because 
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III. LINEAR PERTURBATIONS 
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in this region the massless scalar field is always time-like. To aviod repetition and also for the sake of simplicity, 
we treat only a massless scalar field in this section. The corresponding perturbations for the stiff fluid can be easily 
obtained from the one-to-one correspondence given by Eq. (1.1). With our scope thus focused, we first introduce the 
null coordinates u and v via the relations 

u=-{-u) 2n 1 v = ~(-v) 2n . (3.1) 

Then, we find that the metric in Region II for the background solutions of Eq. (2.2) in terms of u and v takes the 
form ds 2 = 2e 2a dudv - r 2 d6 2 , where 



a(u, v) — 2c 2 ln 



V4 + J8 -l/4 
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r(u, v) = (—u)s(u, v), s(u, v) = 1 — z, z = — > 0, (3-2) 
<f)(u,v) = cln(-u) + ip(z), ip(z) = -2c In (l + z 1/2 ) . (3.3) 

Taking the above self-similar solutions as the background, let us now study their linear perturbations of the form 

F{z,T)=F (z) + eF 1 (z)e kT , (3.4) 

where F = {cr, s, tp}, e is a very small real constant, and r = — ln(— u/uo), with u being constant with the 
dimension of length. Without loss of generality, in the following we set it to 1. Quantities with the subscript "1" 
represent perturbations, and those with "0" represent the self-similar solutions given by Eqs. (3.2) and (3.3). Modes 
with Re(k) > grow with r and are referred to as unstable modes, while those with Re(k) < decay and are referred 
to as stable modes. 

Our treatment of the perturbations is identical to that given in Sec. 4.3 of Rcf. 5), which is referred to as Paper 
I hereafter. The only difference is in the boundary conditions on the SSH, z = 0. As shown in the last section, 
the extension of the background solutions across this surface is analytical only when n is an odd integer. Thus, 
except in the case that n is an odd integer, there is no reason to require that the perturbations be analytical across 
z = 0. Therefore, to impose the boundary conditions properly, we have to distinguish the two cases n = 2m + 1 and 
n ^ 2m+ 1 (to = 0, 1,2, ...). 



Case III. A n = 2m + 1: When n = 2m + 1, the background solutions are analytical across the hypersurface z = 0. 
Therefore, in this case, the boundary conditions are the same as those given by Eqs. (105) and Eq. (108) of Paper I. 
Thus, the analysis given in Sec. 4.3 of Paper I is applicable to this case. In particular, it can be seen that, for any 
given to, these background solutions have N unstable modes, where N is given by JV = n - 3 = 2 (to — 1). Thus, the 
solutions with n = 1 and n = 3 (or m = 0, 1) are stable, while those with n > 5 (or to > 2) are not, and all of them 
have more than one unstable mode. Hence, none of these solutions is critical. 



Case III.B n ^ 2to + 1: When n ^ 2m + 1, the background solutions are not analytical across the hypersurface 
= 0, and there is no reason to require that the perturbations be analytical. Instead, we require that s\, o\ and 
ipi be finite there. In addition, we stipulate that no matter emerge from the self-similar horizon. Thus, we have the 
following conditions: 

si(0), <7i(0), and</?i(0) are all finite, and lim z x ^5£llfj — » o ; (3.5) 

z^o dz 

where \ = c 2 . The boundary conditions at the symmetry axis r = (or z = 1) arc the same as these given in Eq. 
(105) of Paper I. 

Once the boundary conditions are determined, following Paper I we can study the three cases (a) k = 1, (b) 
k = I + 1/2, and (c) k 7^ 1, I + 1/2 separately, where i is a non-negative integer. The perturbations for Case (a) are 
given by Eqs. (110), (111), (114)-(118), while those for Cases (b) and (c) are given by Eqs. (120), (121) and (125) 
of Paper I. Perturbations due to the gauge transformations (79) are given by Eq. (109), which can be obtained from 
Eqs. (120), (121) and (125) by setting ci = = C2, where ci and C2 are free parameters. Thus, to exclude the gauge 
modes from the spectrum of perturbations, in the following we restrict ourselves to the case in which 

|ci| 2 + |c 2 |V0. (3.6) 
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The analyses of Cases (a) and (b) are similar to those given in Paper I. In particular, no perturbations in these two 
cases satisfy the boundary conditions. 

In Case (c), where k ^ 1,1 + 1/2, following the analysis given between Eqs. (120) and (128) of Paper I, one can 
show that the boundary conditions given by Eq. (105) of Paper I at the symmetry axis require c\ = 0, where c\ is 
an integration constant, given in Eq. (125). Then, for s\(z) to be finite at z = it is necessary that Re(k) < 1, as 
one can see from Eq. (120). On the other hand, it can be shown that 
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0, where A(k) = T(k - l/2)/[r(fc)r(l/2)]. Thus, the conditions (3.5) require 
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(3.8) 



Clearly, for any given c (or x = c 2 ), there exist infinitely many modes with Re(k) > that satisfy the conditions (3.6) 
and (3.8). Therefore, all of these solutions are unstable with respect to the linear perturbations. Since each of them 
has an infinite number of unstable modes, they too cannot be critical solutions. 



IV. SUMMARY AND CONCLUDING REMARKS 



In this paper, we have studied a collapsing stiff fluid with HSS in (2 + l)-dimensional circular spacetimes. The 
solutions are usually defined only in Region II [cf. Fig. 1], and Region / is an extended region. The extension is 
analytical only for some of the self-similar solutions. In particular, for the solutions for which n is not an integer, it 
cannot be analytical. 

In the case that n is an even integer, i.e., n — 2m, the extension is analytical, but the stiff fluid has negative energy 
density, and its velocity becomes space-like in the extended region. This is clearly physically unfeasible. Thus, in this 
case, we must remove the requirement of analytical extension. This results in a considerable difference from the case 
of a massless scalar field when we consider boundary conditions on the self-similar horizon for perturbations. Recall 
that, for the scalar field, the extension is analytical, and the only change is that the scalar field becomes space-like 
in the extended region, which is physically feasible. Since the background solutions with n = 2m for a stiff fluid are 
not analytical across z = 0, there is no reason to require that their perturbations be analytical there, in contrast 
to those for the scalar field. Without requiring the analyticity of the perturbations across the SSH, we showed that 
the background solutions with n = 2m for a stiff fluid are unstable and have an infinite number of unstable modes. 
However, in the case of a scalar field, the requirement of analytical perturbations leads to the result that solutions 
with n = 2m have TV unstable modes, where N = n — 3. Thus, the solution with n — 4 has one and only one unstable 
mode, and hence by definition, it is a critical solution for the collapse of the scalar field [5]. Contrastingly, all the 
solutions with either n ^ 2m or n ^ 2m + 1 for both a stiff fluid and a scalar field are non-analytical across z — 0, as 
arc their linear perturbations. As a result, all of these solutions are unstable with respect to linear perturbations and 
have an infinite number of unstable modes. When n = 2m + 1, the extension of the background solutions for a stiff 
fluid is analytical across the SSH. In this case, one may require that their linear perturbations also be analytical. By 
doing so, we found that for any given m, a self-similar solution has N unstable modes, where N = 2(m — 1). Thus, 
the self-similar solutions with m = and m = 1 are stable, while those with m > 2 are unstable, and all have more 
than one unstable mode. 

In review of all the above findings, we can see that all the self-similar solutions, except those with n — 1 and n — 3, 
of a stiff fluid given by Eq. (2.2) are unstable with respect to linear perturbations, and each of them has more than 
one unstable mode. Since a critical solution, by definition, has one and only one unstable mode, we conclude that 
none of these self-similar solutions can be critical for the gravitational collapse of a stiff fluid with circular symmetry 
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in 2 + 1 gravity. This contrasts with the situation for its 3 + 1 dimensional counterpart, where the critical solution 
does have HSS [7,8]. 

Therefore, although between a stiff fluid and a massless scalar field there exists an one-to-one correspondence, 
namely Eq. (1.1), this equivalence is only local, and it is restricted to region(s) in which the scalar field is time-like. 
Globally, the spacetimes for these two cases can be quite different. In the case considered in this paper, this results 
in different boundary conditions, due to their properties in the extended region across the SSH ft is exactly this 
difference that makes the spectra of their perturbations different. 

In the 4D case, similar results were found. In particular, in the collapse of a scalar field, the critical solution is 
periodic and can be time-like only in some regions [1]. Thus, this DSS critical solution is equivalent to a stiff fluid 
only in these particular regions, and it is not a stiff fluid solution in the whole spacetime. For this reason, it cannot 
be critical for the gravitational collapse of a stiff fluid. Instead, it was shown that the critical solution for a stiff fluid 
has HSS [7], subject to the instability with respect to kink perturbations [9]. 

Our results obtained in this paper also clearly show that if critical phenomena arc exhibited in the gravitational 
collapse of a stiff fluid in 2 + 1 gravity, the critical solution must have either discrete self-similarity or no self-similarity 
at all. Certainly, it is also possible that critical collapse does not exist at all in the collapse of a stiff fluid. In any 
case, it would be very interesting to investigate these problems and obtain a definitive answer to the questions. 
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